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^ ■ Abstract 



We have previously found analytically a very unusual and unexpected 
form of confinement in SU(3) Yang-Mills theory This confinement oc- 
curs in the deconfined phase of the theory. The free energy of a single 
static test quark diverges, even though it is contained in deconfined bulk 
phase and there is no QCD string present. This phenomenon occurs in 
cylindrical volumes with a certain choice of spatial boundary conditions. 
We examine numerically an effective model for the Yang-Mills theory 
and, using a cluster algorithm, we observe this unusual confinement. We 
also find a new way to determine the interface tension of domain walls 
separating distinct bulk phases. 
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1 Introduction 



In this paper, we investigate numerically an analytical result discussed in a previous 
paper |1|, where we derived the finite volume dependence of the Polyakov loop in 
cylinders with particular boundary conditions. The analytical result tells us that, 
even in the deconfined phase, an external static quark is still confined. Confinement 
of quarks is usually described by a string tension a, which is the energy cost per 
unit length of the QCD string connecting the fundamental charges. This unusual 
confinement is due to the multiple degenerate deconfined phases and has an effective 
"string tension" a', which depends on the energy cost of a domain wall separating two 
deconfined phases. In this paper, we observe numerically this unusual confinement 
in an effective model for the Yang-Mills theory. Using a cluster algorithm, we can 
accurately determine the energy cost of these domain walls. In addition, we derive 
a new analytical result for the finite volume behavior of the Polyakov loop near the 
deconfinement phase transition, which we also examine numerically with a cluster 
algorithm. In Section 2, we present the analytical calculations, in Section 3, we give 
the details of the numerical work and in Section 4, we briefly draw our conclusions. 



2 Analytical dilute interface gas calculations 

SU(3) Yang-Mills theory at non-zero temperature lives in four-dimensional Eu- 
clidean space-time, with gauge fields A^(x,t) = ieA^x, t)\ a . The non-zero tem- 
perature is T = 1//3, where j3 is the extent in the periodic time direction, i.e. 
A lx (x,t + (3) = A^{x,t). The action S[A^ = J(frft/rf 3 x(l/2e 2 )TrF^F^, is con- 
structed from the field strength F^ u = d^A u — d u A^ + [y4 M ,A„]. The action is 
invariant under gauge transformations A! = g^A^ + <9 M )g, where the SU(3) matri- 
ces g(x,t) are also periodic in the time direction. From our point of view, the most 
important field in the theory is the Polyakov loop, 



$(x) = TrfPexp / dtA 4 (x,t)} 
Jo 



which is complex-valued and is built from the Euclidean time component of the 
gauge field. The Polyakov loop measures the response of the gauge field to the 
presence of a static test quark. Quantitatively, the free energy F of the test quark 
is given by the expectation value ($) = exp(— (3F). If the quark is confined, then 
F diverges and ($) = 0, while if the quark is deconfined, F is finite and ($) ^ 0. 
The Polyakov loop transforms non-trivially under a global symmetry of the action 
S. If the gauge transformations are not quite periodic in the time direction, but are 
instead twisted, i.e. 

g (x, t + (3) = g{x, t)z , ze Z(3) = {exp(27rm/3), n = 1, 2, 3}, (2) 



the action is invariant but the Polyakov loop changes into $'(x) = $(x)z. If ($) = 0, 
then the Polyakov loop value does not break this global symmetry — zero times z is 
still zero. However, if ($) 7^ 0, this global symmetry is spontaneously broken. The 
Polyakov loop is the order parameter for this global symmetry -- the Z(3) center 
symmetry. If the symmetry is intact, the quark is confined; if it is broken, the quark 
is deconfined. In fact, we know that the symmetry is spontaneously broken at high 
temperatures [0. 

We study finite volumes to see whether or not the global symmetry is sponta- 
neously broken in the infinite volume limit. In finite volumes, the system is sensitive 
to the spatial boundary conditions. For example, we could choose all spatial direc- 
tions to obey periodic boundary conditions. However, (<3>) would then always be 
zero, independent of temperature. This is because the Polyakov loop represents a 
static charge and Gauss' law implies that one can't have a charge contained in a torus 
||. We instead apply periodic boundary conditions in only two spatial directions, 
call them x and y. In the third direction z we use what are called charge conjugate 
or C-periodic boundary conditions [Q. If a field is C-periodic, it is replaced by its 
charge conjugate when translated by one period. For example, C-periodicity for 
gauge fields with period L z means that 

A ll (x + L z e z ,t)=A fl (x,t)*, (3) 

where * denotes complex conjugation. Being constructed from the gauge field, the 
Polyakov loop satisfies C-periodicity in the z-direction and we ask that the gauge 
transformations also satisfy this boundary condition. Note that charge conjugation 
is a symmetry of the Yang-Mills action. It is a useful boundary condition to apply 
because the center electric flux coming out of a quark can escape and end on the 
partner anti-quark on the other side of the C-periodic boundary. Now we can have a 
single charge in a finite volume. However, C-periodicity does break the Z(3) center 
symmetry ||, as we see from 

g(x, t)* = g(x, t + (3)*z = g(x + L z e z , t + (3)z = g(x + L z e z , t)z 2 = g(x, t)*z 2 . (4) 

This forces z 2 = 1 and aszG Z(3), therefore z = 1. This means that we no longer 
have the freedom to twist the gauge transformations. This explicit breaking of the 
Z(3) center symmetry vanishes in the infinite volume limit. In a finite volume, ($) 
is always non-zero. If the quark is confined, then (<J>) vanishes in the infinite volume 
limit. Alternatively, if the quark is deconfined, (<3>) remains non-zero, if the infinite 
volume limit is taken in cubic volumes. 

We examine analytically how the Polyakov loop behaves in finite cylindrical 
volumes of dimension L x x L y x L z , where L z 3> L x ,L y . We impose ordinary 
periodicity in the x- and ^/-directions and C-periodicity in the z-direction. First, we 
consider the system at low temperature, where quarks are confined. We view the 
partition function of the system diagramatically as 



exp(-l3f c AL z ), (5) 



where A = L x L y is the cross-sectional area of the cylinder and f c is the temperature- 
dependent bulk free energy density of the confined phase. The partition function is 
the sum of the Boltzmann weights over all possible configurations of the system. At 
low temperature, the entire volume is filled with bulk confined phase. We obtain 
($) diagramatically from 



Z($) = c =exp(-(3f c AL z )Z exp(-f3aL z ). (6) 



Because the Polyakov loop represents a static test quark and the quark is in confined 
bulk phase, there is a QCD string carrying the color flux from the quark to its partner 
anti-quark a distance L z away on the other side of the C-periodic boundary. The 
line in the diagram represents the QCD string, with its string tension o. Knowing 
(<3>) = S exp(— /3aL z ), we find the free energy of the quark is F = — 1//5 In S + crL z . 
This diverges as L z — > oo, i.e. as the quark and anti-quark are pulled infinitely far 
apart — the quark is indeed confined. 

It has been shown that SU(3) Yang-Mills theory has a first order phase transition 
at a temperature T c M. Below T c , external static quarks are confined. Above T c , 
quarks are deconfined. Deconfinement is signalled by a non-zero value for ($) which 
breaks the global Z(3) center symmetry. In the deconfined bulk phase, ($) can take 
any one of three values, 

d>« = ($ , 0), ^ = (-I* , ^$ ), fc® = (-i*o, -^o). (7) 

Each value represents a distinct phase that the system can fall into and these phases 
can be rotated into one another by an element of Z(3). Above T c , the center symme- 
try is broken, giving us three distinct degenerate deconfined bulk phases. Because 
there is more than one phase, different parts of the cylinder can be in different 
deconfined phases. In a cylinder, the phases tend to be separated by interfaces 
which are transverse to the long direction. The phases on either side of an interface 
are labelled by their value of ($), d{ means deconfined phase with the expecta- 
tion value ($) = & l \ An interface between two different deconfined phases has 
a temperature-dependent interface tension, which is the energy cost per unit area 
oidd = F J A. Our diagrammatic expansion of the partition function corresponds to 
summing the Boltzmann weights over all configurations with all possible numbers 
of interfaces, assuming a dilute gas of interfaces @. The diagrammatic expansion 
of the partition function is 



d x 



"2 «3 



d 3 d 2 



+ 



C-periodicity in the ^-direction means that certain configurations are not possible, 
for example, the volume with no interfaces cannot be entirely filled with either 
deconfined phase d 2 or d%. The terms corresponding to the diagrams are 



Z = exp(-f3f d AL z ) 

f Lz 
+ 2 / dz exp(-Pf d Az)j exp(-fia dd A) exp(-f3f d A(L z - z)) + ... 

Jo 
= exp(-pf d AL z )[l + 2 1 exp(-pa dd A)L z + ...]. (9) 

As well as including a Boltzmann weight for the energy cost of the deconfined phase 
with free energy density f d , we must also include the additional energy cost of 
an interface. The factor 7 is included because the interfaces are not rigid and 
actually fluctuate. It has been shown that, in three dimensions, 7 is to leading order 
independent of the area A 0. To sum over all possible configurations, we integrate 
over all possible positions of the interfaces. We sum the interface expansion to all 
orders and find 

Z = ex P (-(3f d AL z + 2 ie xp(-f3a dd A)L z ). (10) 

To calculate (<&) , we again sum over all possible configurations. With the Boltzmann 
weight of each possible configuration, we also include the value of the Polyakov loop 
for that configuration. This gives the expansion 

Z($) = exp(-(3f d AL z ){¥^ + f Lz d Z1 ex V {-(3a dd A) 

Jo 

x —[§W Z + $ i3) (L z -z) + $ {3) z + <$> {2 \L Z - z)\ + ...} 

L z 

= %exp(~Pf d AL z ){l -jexp(-(3a dd A)L z + ...} 

= %exp{-(3f d AL z -^exp{-(3a dd A)L z ), (11) 

where, in the final line, we have summed the expansion to all orders. Dividing out 
the partition function, this gives us the result that 

($) = %exp{-3-fexp{-(3a dd A)L z ). (12) 

This is an unusual result - - it means that the free energy of a static quark is 
F = (—1//?) In $0 + (3'y/f3)exp(—f3a dd A)L z . If A is held fixed, this diverges as 
L z — > 00. As the quark and anti-quark are pulled infinitely far apart, the free energy 
cost becomes infinite. The quark is confined, even though there is no QCD string 
and we only have deconfined phases present. The confinement is due to the multiple 
degenerate deconfined phases lining up along the cylinder in definite domains, which 
are separated by domain walls. We would not see this effect if we had considered a 
cubic volume, as it would have been entirely filled with deconfined phase d\. Also, 
if we had chosen all spatial directions to be C-periodic, the volume would again 
only contain phase d\. The free energy F would then be independent of the volume, 
as we would expect in the deconfined phase. This unexpected confinement has an 
effective "string tension" a' = (37//?) exp(— Pa dd A) which is area and temperature 
dependent. We will observe this unusual confinement numerically and will also use 
the dependence of a' on the area of the cylinder to extract the interface tension a dd . 



We now derive a new analytical result for the finite volume dependence of the 
Polyakov loop near the phase transition. Because the phase transition is of first 
order, as we approach the critical temperature T c from above, the three deconfined 
phases can coexist with the single confined phase. We now implement a different 
choice of spatial boundary conditions. We impose C-periodicity in all spatial di- 
rections. This means that the complex deconfined phases 62 and d% do not appear 
in the cylinder. With these new spatial boundary conditions, the diagrammatic 
expansion for the partition function near criticality is 



d\ 
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(13) 

There are now only interfaces between confined and deconfined phases, which have 
an interface tension a c d and a factor 5 which includes the fluctuations of these domain 
walls. We sum the Boltzmann weights for all possible configurations to obtain 



Z = 2exp(-/3(/ c + f d )AL z /2) cosh(L z \/x 2 + 5- 



(14) 



The quantity x = f3(f c — fd)A/2 measures how far away we are from the finite 
volume critical point (where f c = fd) and 5' = 5exp(— (3a c dA). To calculate Z(§), 
as before, we include the value of the Polyakov loop with the Boltzmann weight for 
each possible configuration. With C-periodicity in all spatial directions, the static 
test quark has a neighboring anti-quark in all directions. If the test quark is placed 
in a region of confined bulk phase, a QCD string connects the quark to the nearest 
anti-quark, which will be in the x- or y-direction. Including the additional energy 
cost of the QCD string, we calculate 



(«> 



{[ 



$n + £' x($ n - £'„* 



exp(L z ^x 2 +5' 2 ) + 



(15) 



x, 



-$o + E' x($ -S' 



cxp( 



-L, 



x/ 



2\x 2 + 5' 



2cosh(L z \/x 2 + 5 ,z ) 



where £[, = S exp(— j3aL±), Lj_ = mm(L x ,L y ). We test this new result in two 
limits. If we are on the deconfined side of the phase transition, then fd < f c and we 
find that, as the volume grows large, (<&) = $o- This is as we expect --a domain 
wall costs more and more energy, so the whole volume is filled with deconfined 
phase d\ and the quark energy is independent of the volume size. If we are on the 
confined side of the phase transition with f d > f c , as the volume becomes large, 
(<&) = S exp(— j3aL±) and F = (— l//3)lnE + aL±. The entire volume is filled 
with confined phase and the quark is connected to the nearest anti-quark by a QCD 
string. This analytical result shows how the finite volume behavior of the Polyakov 
loop changes across the deconfinement phase transition with C-periodicity in all 
directions. We also use this expression to extract the confined-deconfined interface 
tension a c d- 



3 Numerical investigation of the Potts Model 



From the analytical calculations for SU(3) Yang-Mills theory, we expect to see un- 
usual confinement even deep in the deconfined phase. If we can observe such behav- 
ior numerically, we can extract properties of domain walls, such as the deconfined- 
deconfined interface tension add- The analytical calculation assumes that there are 
definite deconfined domains which line up along the cylinder separated by interfaces. 
While this is certainly true if the volume is arbitrarily large, it is not at all clear 
if this phenomenon can be observed in moderately sized volumes. The Polyakov 
loop expectation value should become exponentially small as the cylinder length 
grows. It will certainly be very difficult to detect numerically such a small signal 
for Yang- Mills theory. We want to know if this volume dependence is at all observ- 
able. We will examine the 3-d 3-state Potts model, which is an effective model for 
non-zero temperature SU(3) Yang-Mills theory. We will explain why it is an effec- 
tive model and show how we can observe the unusual confinement even in relatively 
small volumes using highly efficient numerical techniques. 

In the 3-d 3-state Potts model, there are spins <& x at positions x which can take 
one of three values $ x G Z(3) = {exp(27rm/3), n — 1, 2, 3}. The action is 

£[*] = -/*£«•.,•,, (i) 

where the sum is over nearest neighbors. If we globally rotate all spins by <&' x = 
$ x z,z G Z(3), the action is unchanged. This global symmetry is exactly analogous 
to the Z(3) center symmetry of the Yang- Mills theory. Numerical work shows that 



this model has a first order phase transition at a non-zero temperature T c ||12| . For 
T > T c , the system is in a disordered phase, where (<&) is zero and the global Z(3) 
symmetry is intact. For T < T c , there are three ordered phases, where ($) oc 
exp(27rm/3), n = 1,2,3 and the global Z(3) symmetry is broken. The disordered 
phase is exactly analogous to the confined phase of the Yang-Mills theory, while the 
three ordered phases correspond to the three deconfined phases. Because of this 
shared global symmetry and phase structure, the Potts model is an effective theory 
to describe the deconfinement phase transition. We look at the behavior of Potts 
spins in cylindrical volumes which are periodic in x and y and C-periodic in z. We 
can again apply C-periodicity in the ^-direction because charge conjugation is a 
symmetry of the Potts action. Using the dilute gas of interfaces, we calculate that 
the average spin ($) in this volume deep in the ordered (deconfined) regime is 

($) = %exp(-3 1 exp{-pa dd A)L z ). (2) 

Now we have an interface tension add for an interface between two ordered phases 
and 7 accounts for the fluctuation of these non-rigid domain walls. Previously, the 
Polyakov loop fell off exponentially with the cylinder length, signalling confinement. 
This now becomes the exponential fall off of the average spin. Similarly, we calculate 



the average spin near the phase transition applying C-periodicity in all spatial di- 
rections and we see that this depends on a c d, the energy cost of an interface between 
an ordered and a disordered phase. The analytical results are identical, translating 
deconfined and confined bulk phases into ordered and disordered bulk phases. 

We measure the average spin ($) using a cluster algorithm [|IB[]. We pick a spin 
at random from the cylinder as the first element of a cluster. Next, we examine all 
nearest neighbors of this spin. If a neighbor has the same spin value, we include it in 
the cluster with probability P = 1 — exp(— j3). If a neighboring spin has a different 
value, we do not include it - - all the spins in a cluster have the same value. We 
continue to build a cluster until all neighbors of all spins in the cluster have been 
examined. In a multi-cluster algorithm, we keep building clusters until the entire 
volume is decomposed into many clusters. Once all the clusters are built, all the 
spins in each cluster are flipped to a new value. However, it is not always possible 
to flip a cluster. If a cluster winds around a C-periodic direction of the cylinder an 
odd number of times, it must contain spins with real value 1. Otherwise, it would 
not satisfy C-periodicity in this direction. The spins in such a cluster cannot be 
flipped to a new complex value as this would break the C-periodicity. Such a cluster 
is called a wrapping cluster. Any other type of cluster is called non-wrapping and 
can always be flipped. In a multi-cluster algorithm, the average spin is 

< $ > = w( £** ) = ^( E E *. ) = itri Y, N c(m)c ), (3) 

iVy x l\y c x&c l\y c 

i.e. the average spin is the average over all spins in the volume divided by the total 
number of spins Ny. We can sum all the spins in two steps — first we sum over all 
the spins contained in a single cluster C, then we sum over all the clusters. Allowing 
the spins in each cluster to flip to all possible states, we sum the average spin in each 
cluster ((&))c times the number of spins in each cluster Nc- What is the average 
spin in a cluster? In a wrapping cluster, the spins must always have real value 1 
and so these clusters have (($))c = 1- The spins in a non- wrapping cluster can 
take on all possible values and these clusters have (($))c = 0. The entire volume is 
decomposed into clusters, the number of spins contained in the wrapping clusters is 
counted and the clusters are flipped to new values. This is repeated sufficiently often 
to numerically estimate the average number of spins contained in wrapping clusters. 
This is an improved estimator because we only add non-negative numbers in our 
estimate of ($). We can also use a single cluster algorithm. Instead of calculating 
the average spin by summing the average spin of all clusters, we only determine the 
average spin in one single cluster C. The probability of choosing this cluster from 
all of them is proportional to its size, i.e. Pq = Nc/Ny. The average spin is then 

N r ((®))r 

(*> = ( TT V^ > = < « $ »c ) = iVap, (4) 

l\y re 

the sum over all cluster averages is replaced by one cluster average divided by the 
probability Pc of choosing that cluster. The average spin in the entire volume equals 

8 



0.8 



0.6 



0.4 



0.2 



0.0 



-0.2 



-0.4 



-0.6 



-0.8 




0.0 



20.0 



40.0 



60.0 



80.0 



100.0 



Figure 1: The spins Q x are averaged over the transverse plane located at longitudinal 
coordinate z. The plot is of the real and imaginary parts of the planar average of 
the spin versus z. There is one interface, separating two of the ordered bulk phases. 

the average spin in the single cluster that we build, which equals the probability 
-Pwrap that the single cluster is a wrapping cluster. To estimate this probability 
numerically, we build a single cluster; if it is wrapping, we count a 1, if it is non- 
wrapping, we count a 0. We repeat this procedure until we can accurately estimate 
the wrapping probability. Again, this is an improved estimator because we only 
add non- negative numbers in making our estimate of Pwrap- For very long and thin 
cylinders, Pwrap is very small and such a small signal would be very difficult to 
detect if we had positive and negative contributions to our estimator with lots of 
cancellation. The signal is small but the improved estimator does the best possible 
job of numerically estimating the probability. A previous study of the Ising model 
has used an algorithm where boundary conditions can be changed depending on the 



properties of clusters [14 



We first consider the case deep in the ordered phase of the theory, where we 
expect (<&) to decrease exponentially with the cylinder length, which is the signal 
in this effective model for the unusual confinement. Here, we use the single cluster 
algorithm to estimate ($). Previous work has determined that the phase transition 
occurs at (3 C = 0.550565(10) JL2|. The ordered regime is j3 > (3 C . If we go too far away 
from (3 C1 the energy cost of even one interface will be very large, an interface will 
practically never appear in the volume and our analytical calculation will not apply. 
We work in a region (3 C < (3 < 1.005/3 C , where typically there are interfaces in the 
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Figure 2: A plot o/ln($) versus L z , the cylinder length. The lines are the fit to the 
analytical calculation o/ln($) in the ordered phase regime. 



cylinder. The volumes we examine are of transverse dimension L x = L y = 12, 13, 14 
and L z in the range 80 to 120. In Figure 1, we show a typical numerical configuration 
of a volume of dimension 12 x 12 x 100 at a temperature f3 = 1.003/3 C . This particular 
configuration has one interface with bulk ordered phase $ = (—1/2 — iy/3/2)& on 
the left hand side and ordered phase $ = (—1/2 + iy/3/2)$Q on the right. 

In Figure 2, we have plotted the numerical results for ln($) in cylinders of various 
lengths and areas obtained at a particular temperature, here /3 = 1.0035/5 C . We see 
clearly that the average spin falls off exponentially with the length of the cylinder. 
The single cluster algorithm and associated improved estimator are able to detect 
this extremely small signal. We see that our analytical calculation, which is certainly 
true for arbitrarily large volumes, also holds for volumes of moderate size. We fit 
the numerical data to the analytical formula and each line in Figure 2 corresponds 
to cylinders with a fixed area, i.e. 12 x 12, 13 x 13 and 14 x 14. All three lines 
are obtained from one fit, which gives the numerical value of olm{0) with x 2 P er 
degree of freedom ~ 0.7. We take this to be an indication of a good fit. We have 
performed simulations at many values of j3, including one previously examined in 
the literature, allowing us to make a direct comparison of the interface tension. At 
j3 = 0.552, a value of < 
our technique, we find 



0.01796(14) has previously been determined fT5|| . Using 



= 0.01838(42), which is in agreement within error bars. 
In Figure 3, we have plotted the values we have determined for a rfd as a function of 
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Figure 3: The ordered- ordered interface tension add o,s a function of f3. The fit is a 
straight line. 



(3. The data fits very well to a straight line, with x 2 /d.o.f ~ 0.2. The fit is 

a M (P) = -4.54(4) + 8.26(8)0, 



(5) 



where the uncertainties in the final digits are paranthesized. The uncertainties 
appear rather large, but are actually highly correlated. At the critical point, it is 
expected that add(Pc) = 2a c d due to a phenomenon called complete wetting. It is 
tempting to extrapolate our linear fit to the critical point to estimate a c d- However, 
our configurations do not contain disordered phase and our numerical data cannot 
know anything about the energy cost of an ordered-disordered domain wall. We 
cannot extrapolate the fit deep in the ordered phase to the critical point. 

We cannot use the single cluster algorithm to verify our new analytical calculation 
near the phase transition, where we use C-periodicity in all directions. Using the 
single cluster algorithm, ($) is equal to Pwrap- With C-periodicity in all directions, 
the probability of building a wrapping cluster is much greater and so the numerical 
signal is much larger. However, a wrapping cluster cannot be flipped to a new state 
and the configuration is unchanged after this iteration. The more often we build a 
wrapping cluster, the longer it takes for configurations to become independent of one 
another. Thus, a larger signal also means larger autocorrelations. In this case, the 
single cluster algorithm gives much larger integrated autocorrelation times than with 
C-periodicity in the long direction only. This makes the algorithm inefficient and 
we cannot detect finite volume dependence in (<&). In this instance, the multicluster 
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Figure 4: A plot of ($) versus (3 in cylinders of length L z = 80 and areas from 
12 x 12 to 16 x 16. The curves are the fit to the analytical behavior of ($) across 
the phase transition. 

algorithm is far superior. By decomposing the entire volume into clusters, there 
will always be many non-wrapping clusters which can be flipped and configurations 
very quickly become decorrelated. Using the multicluster algorithm, we examine 
cylinders with dimensions L x = L y = 12, 13, 14, 15, 16 and L z = 60, 70, 80, 90 and at 
temperatures in the range (3 G [0.545, 0.551]. We fit all of the data to the predicted 
analytical behavior of ($). In Figure 4, we plot some of the numerical data and 
the fit of ($) as a function of (3. From the fit of all the data, we find the ordered- 
disordered interface tension is a c d = 0.0015(2), with % 2 /d.o.f ~ 1.1. Translated into 
our notation, a value of a C( i = 0.00148(2) has previously been found [0. From the 
fit, we also obtain j3 c = 0.55070(17), in comparison to the previously quoted value 
(3 C = 0.550565(10) [0. These values are consistent and we feel that the relatively 
large percentage errors in our estimates of the interface tension and the critical 
temperature may be because we work in smaller volumes. There are also other 
techniques to extract the interface tension near criticality [p 



4 Conclusions 



Analytically, we predicted that a very unusual type of confinement occurs even in 
the deconfined phase of SU(3) Yang-Mills theory. The free energy of a single static 

12 



quark diverges, even though it is contained in deconfined bulk phase. The signal for 
this unusual confinement is that the expectation value of the Polyakov loop becomes 
exponentially small. We have successfully observed this phenomenon numerically 
in an effective theory. It was not a priori clear that the exponentially small signal 
could be measured or that this confinement would be observable in moderately sized 
volumes. The single cluster algorithm and associated improved estimator can detect 
this very small signal in reasonably sized cylinders and accurately determine the 
energy cost of domain walls between distinct bulk phases. Our new analytical result 
predicts the finite volume behavior of the Polyakov loop across the deconfinement 
phase transition and we observe this behavior numerically, allowing us to estimate 
the energy cost of domain walls at criticality. Interestingly, this is an example where 
a multicluster algorithm is far superior to a single cluster algorithm. 



5 Acknowledgements 

I would like to acknowledge the essential contribution of U.-J. Wiese to this project. 

References 

[1] K. Holland and U.-J. Wiese, Phys. Lett. B415 (1997) 179. 

[2] L. D. McLerran and B. Svetitsky, Phys. Rev. D24 (1981) 450. 

[3] E. Hilf and L. Polley, Phys. Lett. B131(1983) 412. 

[4] A. S. Kronfeld and U.-J. Wiese, Nucl. Phys. B357 (1991) 521. 

[5] U.-J. Wiese, Nucl. Phys. B375 (1992) 45. 

[6] R. V. Gavai, F. Karsch and B. Petersson, Nucl. Phys. B322 (1989) 738; 
M. Fukugita, M. Okawa and U. Ukawa, Phys. Rev. Lett. 63 (1989) 1768; 
N. A. Alves, B. A. Berg and S. Sanielevici, Phys. Rev. Lett. 64 (1990) 3107. 

[7] B. Grossmann, M. L. Laursen, T. Trappenberg and U.-J. Wiese, Nucl. Phys. 
B396 (1993) 584. 

[8] V. Privman and M. E. Fisher, J. Stat. Phys. 33 (1983) 385; 

E. Brezin and J. Zinn- Justin, Nucl. Phys. B257 [FS14] (1985) 867. 

[9] Z. Frei and A. Patkos, Phys. Lett. B229 (1989) 102. 

[10] T. Trappenberg and U.-J. Wiese, Nucl. Phys. B372 (1992) 703. 



13 



[11] K. Kajantie, L. Karkkainen and K. Rummukainen, Nucl. Phys. B333 (1990) 
100; Nucl. Phys. B357 (1991) 693. 

[12] W. Janke and R. Villanova, Nucl. Phys. B489 [FS] (1997) 679, and references 
therein. 

[13] R. H. Swendsen and J. S. Wang, Phys. Rev. Lett. 58(1987) 86; 
U. Wolff, Phys. Lett. B228 (1989) 379. 

[14] M. Caselle and M. Hasenbusch, Nucl. Phys. B470 (1996) 435. 

[15] M. Caselle, F. Gliozzi, P. Provero and S. Vinti, Nucl. Phys. Proc. Suppl. 34 
(1994) 720. 

[16] M. Schmidt, Z. Phys. B95 (1994) 327; 

J.-D. Wang and C. DeTar, Phys. Rev. D47 (1993) 4091; 

B. A. Berg and T. Neuhaus, Phys. Rev. Lett. 68 (1992) 9; 

B. Grossmann, M. L. Laursen, T. Trappenberg and U.-J. Wiese, Phys. Lett. 

B293 (1992) 175. 



14 



